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Abstract 

We study the large-time behaviour of the solutions u of the evolution equation involving nonlinear diffusion 
and gradient absorption 

afM-ApW+|Vu|9 = 0. 

We consider the problem posed for x € and t > Q with non-negative and compactly supported initial 
data. We take the exponent p > 2 which corresponds to slow p-Laplacian diffusion, and the exponent q 
in the super linear range 1 < g < p — 1. In this range the influence of the Hamilton- Jacobi term iVul*? is 
determinant, and gives rise to the phenomenon of localization. The large time behaviour is described in terms 
of a suitable self-similar solution that solves a Hamilton- Jacobi equation. The shape of the corresponding 
spatial pattern is rather conical instead of bell-shaped or parabolic. 

1 Introduction 

Researchers have been interested for decades in the long time description of evolution processes 
where diffusion is combined with other effects, notably reaction, absorption and/or convection. The 
equations under study are evolution equations of parabolic type, mostly nonlinear and possibly 
degenerate parabolic. A general form of such equations is 

N 

dtu = '^dx^Ai{t,x,u,Vu) + B{t,x,u,Vu), (1) 

i=l 

where the Ai and B are nonlinear functions with suitable structure conditions as in j30j . and Vu 
denotes the spatial gradient of u. Under such a generality one cannot be very specific about the 
concrete asymptotic behaviour of the solutions as t ^ oo. On the contrary, when one concentrates 
into specific forms for Ai and B, and also imposes specific boundary conditions, or works in the 
whole space for concrete classes of initial data, then much is known about the long time behaviour 
of the solutions, and the different types of possible behaviours have been carefully classified, or are 
in the process of being classified. 
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Following this spirit, we will concentrate here on equations where the right-hand side of the 
equation ([1]) combines nonlinear diffusion of the p-Laplacian type with nonlinear absorption of 
Hamilton- Jacobi type. We pose the evolution on M^, > 1, as spatial domain, and take as initial 
condition a non-negative and integrable function uq decaying sufficiently rapidly to zero as \x\ oo 
(we may assume uq to be compactly supported and bounded to make the analysis simpler). Under 
these assumptions, we want to classify the different types of long type behaviours that may arise 
from the interaction between the two physical effects at work. As a new contribution to the topic, 
we describe in detail a situation in which the evolution leads to localized patterns with a precise 
power-like decay in time and conical spatial shape. 

1.1 The case of linear diffusion 

Before we enter into that study and in order to motivate the issue, it will be useful to briefly 
recall the main results in the case of standard linear diffusion. The large time behaviour of non- 
negative and integrable solutions to the diffusive Hamilton-Jacobi equation (also called viscous 
Hamilton-Jacobi equation) 



has been studied recently by several authors. It turns out that such behaviour is not unique and 
strongly depends on the value of the parameter q € (0,oo). More precisely, always for initial data 
decaying sufficiently rapidly to zero as — > oo, the analysis reveals the existence of two critical 
exponents qi := 1 and (7* := (A^ -|- 2)/{N + 1), and three different types of large time behaviour 
corresponding to the parameter intervals q € (0, gi), q G {qi,q*), and q € (g*,(X)), respectively. To 
this we must add the study of the critical cases q = 1 and q = q*. 

To be more precise, in the upper range q > q^, the large time dynamics is governed by the sole 
diffusion term and u behaves as a multiple of the fundamental solution to the linear heat equation 
im I17|. After convenient renormalization we thus see the typical bell-shaped Gaussian proffie, in 
the sense that 



what is called in [M] asymptotic simplification, and it also happens for the semilinear equation 
dtu - Au -u"" = in Q when r > := {N + 2)/N. 

In the intermediate range q G (1, g*), the large time behaviour results from the combined effects of 
the diffusion and absorption terms and is fully described by the self-similar Very Singular Solution 
to equation (j2]), cf. [H]. The pattern that we see after the corresponding renormalization, Fq{y), 
is a modification of the Gaussian profile corresponding to the Very Singular Solution and changes 
with q. 

As usual in dynamical studies, the large time behaviour for the critical exponents qi and is 
peculiar. For g = it is still given by a multiple of the fundamental solution to the linear heat 
equation but with extra logarithmic factors resulting from the fact that the dynamics already feels 
the effects of the absorption term [22]. This is usually referred to as resonance. 

On the other hand, when q equals qi = 1 then ([2]) is still a nonlinear equation but with the same 
homogeneity as a linear equation, and the large time behaviour does not seem to be thoroughly 





with y = x/t 
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understood. Nevertheless, there are several possible temporal decay rates for u depending on the 
initial condition [13\ I14j. This is an important issue that needs further investigation. 

Finally, when the parameter q is below qi the nonlinear absorption term becomes dominant and 
diffusion plays a secondary role for large times. The typical phenomenon of the lower range is finite 
time extinction, which takes place for q G (0,gi) |12l I23| . Summarizing, the nonlinear absorption 
term only prevails for large times if (7 € (0, gi). But the exponent = 1 is actually somehow 
"doubly" critical: indeed, not only does the nonlinearity rule the dynamics but it is also no longer 
locally Lipschitz continuous for q G (0,gi), thus giving rise to singular phenomena such as finite 
time extinction. 



1.2 Nonlinear diffusion. Localized non-diffusive patterns 

In an attempt to elucidate the true role of the absorption term when it governs the dynamics, we 
devote this paper to investigate in detail the large time behaviour of non-negative solutions u to 
the related Cauchy problem 

dtu- ApU + \Vu\'^ = 0, {t,x)eQ, (3) 
ti(0) = no > , X G , (4) 

where the linear diffusion operator Au is replaced by the p-Laplacian operator 

Apu := div {\Vu\p~^ Vu) . (5) 

When p > 2, ^ is a. quasilinear degenerate parabolic equation which reduces to the semilinear 
equation ([2]) for p = 2, and has been studied recently: the existence and properties of self-similar 
solutions are investigated in [251 E3], while qualitative properties of non-negative and integrable 
solutions are obtained in [H \TU[ I37j . In particular, two critical exponents have been recently 
identified for ([3]) [U [lO] , namely 

N 

qi := p — 1 and := p 



N + 1 



As for p = 2, the nonlinear absorption term governs the large time dynamics for q G (0, qi) but we 
now have qi > 1 since p > 2. However, the singular phenomena like extinction in finite time are 
only expected to happen in a still lower subregion q G (0, 1). 

Consequently, in the new intermediate range q G (1, 9i) that we want to explore here the nonlin- 
earity is locally Lipschitz continuous and no singular phenomenon (such as extinction) can occur. 
Still, the absorption term strongly affects the evolution for q G {l,qi). Two results are established 
in |10j in that parameter interval for compactly supported initial data: there are positive constants 
C > and R> (possibly depending on the initial data) such that 

Vu{t) C B{0,R) and ||n(t)||oo < C t-Vte-i) 

for alH > [here and below we use the notation Vuit) = {x G : u{t, x) > 0| for the positivity 
set of the function u at time t and we use the abbreviated notation u{t) instead of u{t, •) and so 
on, whenever there is no fear of confusion]. 
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These properties are reminiscent of those enjoyed by the viscosity solutions to the Hamilton-Jacobi 
equation 

dth+\Vh\'' = in Q (6) 

with compactly supported initial data. Indeed, denoting by h the viscosity solution to equation ^ 
with initial condition h{0) = u{0), we have for each t >0 (cf. [7]) 

sup|ii/(«"^M|/i(i)lloo| <oo, and Vh{t) = Vh{0). 

This is in sharp contrast with the behaviour of the solution w to the pure p-Laplacian equation 
dtw — ApW = in (0, oo) x with the same initial condition w{0) = u{0). Indeed, the L°°-norm 
of w{t) only decays at the much slower rate 

||u;(t)||oo = Or^/(^(^-')+^)), 

and the positivity set Vw{t) = {x E : w{t,x) > 0} of expands and as t ^ oo it fills the 
whole space M.^ . The absorption term |Vti|'' in ([3]) thus prevents this expansion if q G and 
the positivity set of u{t) remains localized in a fixed ball of for all times. 

The property of localization is already known to be satisfied by compactly supported non-negative 
solutions to second-order degenerate parabolic equations with a sufficiently strong absorption in- 
volving the solution only as, for instance, dtW — ApW + = in Q when r G — 1), cf. 
[211 EH [37] . It had apparently remained unnoticed for second-order degenerate parabolic equations 
with an absorption term depending solely on the gradient. 

According to the previous discussion, the qualitative results obtained so far on compactly sup- 
ported non-negative solutions to ([3]) show evidence of the domination of the nonlinear absorption 
term for large times. The purpose of this paper is to go one step further in that direction and show 
that u{t) behaves like a self-similar solution to ([6|) as t — > oo. 

For comparison, it is interesting to notice that this localized asymptotics is quite different in shape 
from the behaviour of the diffusion-absorption equation dfW — Aw"^ -|- = in Q that offers a 
strong similarity in other respects with our present problem when 1 < r < m . The corresponding 
localized behaviour has been studied in [15^ [TBI [H] ^ind leads to a mesa-like pattern. In this paper 
we show that a conical pattern Voo is gradually formed, precisely given by formula (jl5p . and is 
actually a viscosity solution to the stationary equation 

\S/v\'^-v = in . (7) 

It is associated to an exact self-similar viscosity solution of the Hamilton-Jacobi equation ([6]). 

As a final contribution, we investigate conditions under which the support does not move at all, 
so-called inGnite waiting time, so that what happens for the whole time span is a process of internal 
reorganization. We refer to [1] and [36] for waiting times in porous medium flows. In those cases 
the waiting time is always finite. An infinite waiting time is described in [181 119] for the equation 
dtw — Aw^ + = in Q when 1 < r < m. Properties of localization and waiting times for more 
general equations of the form ([T|) are considered for instance in [21 [31 [21] . 
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2 Preliminaries and main results 



Before stating our results, let us first specify our assumptions and recall the properties of solutions 
to dS]), ([3]) with non-negative and compactly supported initial data established in [10]. We assume 
that 

p>2 and qe{l,p-l), (8) 
and that the initial condition uq enjoys the following properties: 

uo£W^'°°{R^), uo>0, Vo:= {x eR^ : uo{x) > 0} C B{0, Rq) , uq ^ , (9) 

for some Rq > 0. Let us also introduce an exponent that will play a role in what follows: ^ := 
l/(g(iV + 1) - N) > 0. It goes from ^ = 1 for g = 1 to ^ = l/{N{p - 2) + p - 1) when q = p - 1. 

Proposition 1 Under the above assumptions, the Cauchy problem ([3]), ^ has a unique non- 
negative viscosity solution 

u e i3C([0, oo) X M^) n L°°(0, oo; VF1'~(M^)) 

which satisfies: 

(i) There is Roo > such that 

supp {u{t)) C B{0,Roo) for all t>0. (10) 
(a) There are positive constants Cq and Ci such that 

" Mt)\\oo<Co Ms)\\f {t-sr^i, 

(11) 

ll^^Wlli < Ikolli , ||'"(i)lloo < Ikolloo , ||Vu(t)||oo < llVuolloo , (12) 
(\\u{t)\\i + \\u{t)\\^ + \\Vu{t)\U < Ci (13) 

for allt>0 and s G [0,T). 



Here and below BC{[0,oo) x M^) denotes the space of bounded and continuous functions on 
[0, oo) X and || • \\r denotes the L''(M^)-norm for r G [l,C)o]. 

Assertions (fTO]l . (fTT]l . and (fT2l) are proved in [ini Theorem 1.5 (i)] and Proposition 1.3], 
respectively, while (I13p is established in [101 Corollary 1.6 (i)] for the L^-norm. The estimates ()13p 
for ||ti(t)||oo and ||Vti(t)||oo follow then from that for ||ti(t)||i by (fTT|) with s = t/2. 
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Profile V for q=5/4 



ProfileVforq=10 





Figure 1: Profile Voo for q = 5/4 (left) and q = 10 (right) 
2.1 Statement of the main asymptotic results 

As indicated above we call Vuit) the positivity set of the solution u at time t > 0, that we abbreviate 
as V{t). Since u{t) is compactly supported and continuous for each t > by Proposition [U 'P(t) is 
a bounded open subset of for all t > 0. Much more can actually be said on the family {V{t)}t>o 
and it will allow us to identify the large time behaviour of u. 



Theorem 2 The mapping t i — > V{t) gives a non- decreasing family of bounded open subsets q 
and moreover 

Voo '■= ^(^) ^-5 bounded open subset of . 



In addition. 



where 



t>o 



lim 

t— >oo 



0-1 / \ 9/(9-1) 

^oo(x):=^ ynf^JI.-l}) for x . 



(14) 
(15) 



We first point out that U{t,x) := t-V(<?-i) (x) is a self-similar (viscosity) solution to ^ which 
vanishes outside Voo and is positive in Vcxi- Therefore, p4|) asserts that the diffusion term no 
longer matters for large times as claimed previously. Still, it might influence the shape of Voo 
during the time evolution. We will come back to this issue in Theorems [3] and H] below. Also, 
the final (rescaled) pattern Voo has a conical shape as illustrated in Figure 1 for N = 2 and 
Voo = {{x,y) G [-2,3] X [0,9] : ^2 < y < x + 6}. 

The proof of Theorem [2] is divided in several steps and is performed in the forthcoming Sections [3] 
and m We nevertheless outline its main steps now: we first establish the monotonicity property of 
the family {'P(i)}t>o by comparison arguments, and lower bounds for u{t, x) as well. We next turn 
to the proof of (|14p and introduce the self-similar variables; we keep the space variable x (since we 
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expect localized behaviour), and (as usual) introduce logarithmic time 



r :=ln(l + (g-l)t)/(g-l), 

as well as the new unknown function v = v{t, x) defined by 

u{t,x) = {l + {q-l)ty^/^'^^^^ v{t,x) , (i, x) € [0, oo) X . (16) 

Usually, v{t) = v{t, •) is called the renormalized trajectory. Clearly, v solves 

drV + \Vv\'^-v = Apv , {T,x)eQ, (17) 

v{0) = uo, x€M^. (18) 

The convergence (fT4l) then reads \\v{t) — Voo||oo — > as r — > co and we are thus left to establish 
the convergence of u(r) as r — > oo. We first notice that, as q < p — 1, the coefficient in front of 
the diffusion term vanishes for large times, so that the diffusion term is clearly negligible for large 
times. To identify the large time behaviour of v{t) we use the half-relaxed limits technique [S] in 
the same fashion as in [SH Theorem 1] and [32l Section 3], and introduce 

v^{t,x) := liminf v{a + s,y) and v*{t,x):= limsup v{a + s,y) 

for (r, x) € Q. We first prove that f * and v* do not depend on r and are viscosity supersolution 
and subsolution, respectively, to the Hamilton- Jacobi equation |Vz|'^ — z = in M^. Noticing that 
and V* are positive on the same open subset -Poo of M^, a comparison argument can be applied 
to conclude that v^, = v*, from which the convergence of f (r) follows. 

2.2 The question of waiting times 

In Theorem [2] we have identified the first term of the expansion of u{t) as t ^ oo. Still there is 
an unspecified quantity, namely the set Voo on which additional information would be welcome. In 
particular, is there a class of initial data for which the positivity set does not change during time 
evolution? The next result gives a positive answer to this question. 

Theorem 3 Assume that there exists 6 > such that 

uo{x) < ao |x - xol^^"")/^^"^^"^ for x G B{xo, 6) and xq G dVo , (19) 
llnolloo <ao ^(^-''^/(P-''-^), (20) 

where 

ao:=^ ^ iV-l + ^- . (21) 

p-q V p-^-qJ 

Then V{t) = for every i > 0, and thus Voo = "Po- 
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When a point of the boundary of the positivity set stays as such for a time we say that there is a 
waiting time at that point. Waiting times are typical in degenerate parabohc equations with slow 
diffusion, like dtw = Aim™ for m > 1 or dtw = ApW for p > 2, cf. [36], but in these cases they are 
finite. In the present case we have exhibited an infinite waiting time at all points of the boundary 
of the initial support. 

For small values of q, the behaviour reported in Theorem [3] ceases to be observed as soon as ()19p 
is not fulfilled, showing the optimality of the exponent {p — q)/ip — ^ — q) for the support not to 
evolve. 

Theorem 4 Assume that xq £ dVo and there is 6 > such that 

|(p-g)/(p-i-g) x£B{xo,6). (22) 

iV + 1'2, 

there is Aq > depending only on N, p, and q such that, if A > Aq, 

u{t, xq) > for all t > 0. 
In particular, xq £ 'P{t) for each t > 0. 

A similar result is actually also valid for q £ [q2,P — ^) but requires a stronger assumption than 
(j22p . see Proposition [TT] below. 



q < q2 :=mmip- ——,-)■ (23) 



3 Time monotonicity of the positivity set 

We first establish the time monotonicity of V{t). 

Proposition 5 For ti G [0, oo) and t2 S (ti,oo) we have 'P(ti) Q "^(^2) ^.i^-d 

Voo '■= 'Pit) is a bounded open subset of . (24) 
t>o 

The proof relies on the construction of a suitable subsolution which we perform next. 

Lemma 6 Fix R> and put 

SA{t) ■.= A{1 + t)-V{'/-i) (ij2 _ |2.|2)9/{g-i) ^ ^ ^ ^^Q^ ^ 

There exists Ar > such that sa is a (viscosity) subsolution to ([3]) in (0, 00) x B(0,R) for every 
Ae[0,AR]. 
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Proof. Introducing a{x) := {R^ - jxp)^/^^ for x S B{0,R), we observe that a and |Vc7|p ^Vct 
both belong to C^{B{0,R)). Therefore, if {t,x) G [0, oo) x B{0,R), 

^ ^ {dtSAit, x) - ApSAit, x) + \VsA{t, x) 



^ cj(x) 

+ (l+t)(p-l-g)/(</-l) ^^"^^ 1^1 



. 1 \q-i 

+ ' (iV + P - 2) ^{2{p-l-g)/{c-l))+p-2 < Q 

if € [0, Aii\ where is sufficiently small and depends only on A^, and R. □ 

Proof of Proposition [5], Fix ti > and xi G V{ti). Owing to the continuity of a; i-^ u(t, x) there 
are > and > such that u{ti,x) > 5 > for x G B{xi,R). Choosing A G {0,Ar) such that 
A<5 i2-29/(9-i)^ we notice that 

SA{ti,x-xi) <A R^^/^'i^^^ <5 <u{ti,x) for xeB{xi,R), 

and 

X — xi) = < x) for (t, x) G [ti, oo) X 9i?(xi, i?) , 

the parameter A/j and the function sa being defined in Lemma[6l By Lemma[6l {t, x) i — > SAit, x — 
xi) is a subsolution to ^ in [ti,oo) x B{0,R) and we infer from the comparison principle [201 
Theorem 8.2] that SA{t,x — xi) < u{t,x) for (t, x) G [ti,oo) x B{0,R). In particular, n(t, xi) > 
SA(t,0) = A (1 + t)-V(9-i) ^2q/(q-i) > Q for t > ti. Therefore, xi G for every t > h, which 
proves the first assertion of Proposition [5j The set "Poo is then clearly an open subset of which 
is bounded hy □ 

Another useful consequence of Lemma [6] is a lower bound for u{t) which implies in particular that 
|[u(t)||oo cannot decay to zero at a faster rate than for large times. 

Corollary 7 For each x G Voo there are Tx > and > such that 

u{t, x) + for t>Tx. (25) 
In particular, there exists C2 > such that 

h(i)||oo >C2 t-^/^"-^) for t>l. (26) 

Proof. Consider x G Voo- Then there is T^; > such that x G V{t) for t > T^. The continuity of 
y I — > u{Tx, y) ensures that there are 6 > and R > such that u{Tx, y) > (5 > for y G -B(.x, i?). 
Choosing A G (0,Ajj) such that A < 6 R"'^i/(i~^) we argue by comparison as in the proof of 
Proposition [5] to conclude that SA(t,y — x) < u{t,y) for {t,y) G [^^,00) x B{x,R) (recall that Ar 
and SA are defined in LemmaE]). Consequently, u(t,x) > SA{t,0) = A (1 + t)-^/(^~i) /?2g/{g-i) > q 
for t > Tx, whence (j25p with Sx = A R'^'^^^'^~^\ The lower bound (j26p is then a straightforward 
consequence of (f25|) with any x G "Pq- D 
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4 Convergence to self-similarity 



We now investigate the properties and convergence for large times of the function v defined in (jl6p . 
For r > let Vv{t) := {x G : v{t,x) > O} be the positivity set of v{t). According to (fT6|) . 
Vvir) is given by 

/e(9-i)r _ i\ 
V.v(t) =Vu{ ^ _ ^ J for r>0. 

We next state some bounds for the rescaled function v which are easy consequences of (jlOp . ()12p . 
(dl, and ([25]). 

Lemma 8 There is C3 > such that 

lb(r)||oo + ||V^;(r)||oo <C3 for r>0. (27) 
In addition, for each x € 'Poo; there are Tx and > such that 

v{t) >ex>0 for t > . (28) 

Finally, r 1 — > VviT) is a non- decreasing family of bounded open subsets o/M^ and 

U rv{T)= Voo, 
r>0 

the limit set Voo being defined by (124p . 

Proof. The estimates (I27p readily follow from (I12p . (jl3p . and (jlGh . while ()28p and the properties 
of {'Pi;(t)}t->o are straightforward consequences of Proposition [5] and Corollary [71 □ 

We turn to our main task, i. e., the behaviour of v{t) as r — > 00, and actually aim at showing the 
convergence of v{t) as r — > cxd to an asymptotic profile. For that purpose we use the half-relaxed 
limits technique [9] in the same fashion as in [321 Section 3]. For (r, x) E [0, cxo) X we define 

v^,{x) := liminf v{a + s,y) and v*{x):= limsup v{a + s,y) (29) 

(fT,j/,s)^(r,x,oo) ((T,j/,s)^(r,x,oo) 

and first note that the right-hand sides of the above definitions indeed do not depend on r > 0. In 
addition, 

Q <v^{x) <v*{x) for xGM^ (30) 

by ()29p . while Lemma [8] and the Rademacher theorem clearly ensure that and v* both belong 
to iyi'°°(]R^). Finally, by [8l Theoreme 4.1] applied to equation ([IT]), v* and v^, are viscosity 
subsolution and supersolution, respectively, to the Hamilton-Jacobi equation 

i?(z, Vz) := iVzl" - z = in . (31) 
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The next step is to show that v* and actuahy coincide. At this point we emphasize that v* 
and enjoy an additional property, namely, 



v*{x) = v^{x) =0 for X G \ 
v*{x) > > for X & Voo 



(32) 



(a similar situation is encountered in the proof of |3H Theorem 1]). Indeed, if x "Pooi the time 
monotonicity of r Vv{t) (Lemma [S|) warrants that x ^ Vv{t) for all r > 0. Then, if s„ — > oo, 
an 0, and x„ ^ x are such that u((7„ + s„, x„) — > v*{x) as n ^ cx), it follows from (p7|) that 



■y(0"„ + Sn, < ^((Tn + S„, x) + |x - X„| || V7;(o-„ + Sn) 



< Co |X - Xr. 



Passing to the limit and recalling ([50]) gives the first assertion in ([52]) . Consider next x G "Pqo and 
recall that there are and e^; > such that v{t,x) > Ex for r > by (pHj) . Pick sequences 
(o-n)n>i, (sn)n>i, and (x„)„>i such that cr„ ^ 0, s„ ^ oo, x„ ^ X, and i;(<t„ + s„,x„) — > f*(x) 
as n — > oo. For n large enough we have o"„ + s„ > and we infer from Lemma [8] that 

V{(7n + Sn,Xn) > v{an + Sn,x) - C3 |x - Xn\ > Er, - C3 \x - X„| . 



Letting n — > oo gives u*(x) > Sx > and completes the proof of (|32 
We next introduce 

Q 



V,{x) :-- 



q-1 



,(x)(''-^)/« and V*{x) :-- 



1 



v*{x) 



(33) 



for X G 'Poo- Arguing as in [21 Corollaire 2.1] or [BJ Proposition II.2.5], it easily follows from the 
properties of v* and that V* and are viscosity subsolution and supersolution, respectively, to 
the eikonal equation 



\Vz\ -1 = in Vo 



(34) 



Now, Hei depends solely on Vz and is a convex function of Vz, so that the assumptions (HI), 
(H2), and (H4) in [26] are clearly fulfilled. Furthermore, Hei{0) = — 1 < which warrants that 
the assumption (H3) in [26] is also fulfilled (the function ip in [26] being here identically zero). 
Since K = ^* = on dVoo, we are in a position to apply [261 Theorem 1] to conclude that 
V*{x) < 14 (x) for X G Poo- Recalling ([32]) and §2l) we end up with V*{x) = V^{x) for x G Poo and 
thus v^:{x) = v*{x) in M^. Owing to [H Lemme 4.1] or [6l Lemma V.1.9], the equality = v* and 
(j29p provide the convergence of {v{s))s>i towards uniformly on every compact subset of as 
s ^ oo. Since all these functions are compactly supported in the closure of Poo, we conclude that 



lim \\v{t) — 

T— ►oo 

Returning to the original variables (t, x), we obtain 



0. 



whence 



lim 

t— >oo 



lim 

t^oo 



il + iq-l)ty/^'^~^\{t)-v 

tl/fa-l)^(t)_(g_l)-l/(5-l) 



0, 



0. 



(35) 
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It remains to identify v^, or equivalently V*. Since the latter is a viscosity solution to |VKi,| — 1 = 
in Voo with 14 = on dVoo by the previous analysis, we infer from [6l Remark II. 5. 10] that 

V^(x) = inf {\y — x\} = inf {ly — a:!} for xGVoo- 
Consequently, = Voo, the latter being defined by (fT5]) . and (fT^ follows from ([35]) . □ 

5 Invariant and moving posit ivity sets 

We first show that, if uq vanishes sufficiently rapidly near a point xq, then u{t, xq) = for all t > 0. 

Lemma 9 Assume that there are xq € and 6 > such that 

|(p-g)/(p-i-g) x£B{xo,6), (36) 

\\uo\\oo<ao6^P~''y^P-'^~'\ (37) 
the parameter ao being defined by (pT|) . Then, 

u{t,xo) = for t>0. (38) 

Proof. We adapt an argument from \29\ Theorem 8.2] and put Si{x) := ao \x — xo\^p~'^'^^^p~^~'^^ 
for x G M^. An easy computation shows that 

-Ap5i(x) + |V5i(x)|« = for x gR^ . 

Furthermore, we have uo{x) < Si{x) for x £ B{xo,6) by ([36]) and u{t,x) < ||mo||oo < Si{x) for 
{t,x) G [0, oo) X dB{xo,6) by 0^ and (f37|) . The comparison principle [201 Theorem 8.2] then 
entails that u{t,x) < Si{x) for {t,x) G [0, oo) x i?(j;o,(5). Consequently, u{t,xo) < Si{xo) = for 
t > and the lemma is proved. □ 

Proof of Theorem [S], We first consider xq G OVq. Owing to (fT9]) and ([20]) we are in a position 
to apply Lemma [9] to deduce that u{t, xq) = for all t > 0. We next consider xq lying outside 
the closure of Vq and x G i?(xo,(^). Either x ^ Vq and tiol^;) = < Cq ja; — xqI^^"''/^^"^"''^. 
Or X G 7^0 and there exists -& G (0, 1) such that x.^ := -dx + {1 — '&)xo belongs to dVo- Then 
\x — x^\ = (1 — ??) \x — xo\ < 6 and we infer from (119p that uq{x) < qq \x — x^\^^~'^^^^~^~'^^ < 
ao \x — xo\^~'^^^P^^~^'\ The assumptions of Lemma[9]are then fulfilled by uq at xq and we conclude 
again that u{t, xq) = for all t > 0. 

We have thus shown that V{t) C Vq for t > 0. The opposite inclusion readily follows from the 
time monotonicity of the positivity set established in Proposition [5j □ 

We now turn to the proof of Theorem H] and first establish the following result, which is in the 
spirit of [5l Proposition 3.2]. 
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Lemma 10 Consider q G (1,^2), q2 being defined in (|23|) . There are positive real numbers C4 and 
C5 depending only on N , p, and q such that, if uq is an initial condition fulfilling ([9]) and 



M := uoix) dx>C5, 

Jb{o,i) 

then the corresponding solution u to ^ satisfies 



n(l,0) > Ca AF'' with T] : ^ 



N{p -2)+p 



Proof. We adapt the proof given in [35] for the porous medium equation, the main difference 
being that the L^-norm of u is not constant in our case. In the following, we denote by Cj, i > 6, 
positive constants depending only on A'^, p, and q. 

Let us first assume that supp (uq) C B{0, 1), so that ||tio||i = M. Denoting by w the solution 
to the p-Laplacian equation dtw — ApW = in Q with initial condition w{0) = uq, the comparison 
principle entails that u{t,x) < w{t,x) for {t,x) € [0, 00) x M^. By [28^ Proposition 2.2] and [W\ 
Proposition 1.3] we have 

II^WIloo < |b(t)||oo<C6M^"'t-^\ 

u{t, x) = v{t, x)=0 if |x| > 1 + Cr M'^^P'^^ t"^ , 
||Vn(t)||oo < C M^" r^^^+i) , 

for t > 0, the constants Cq and C7 depending only on N, p, and q. We next use the reflection 
argument of Aleksandrov as in the proof of [5l Lemma 2.2] to deduce from ([3]) that 

u{t,0) >u{t,x) for t>0 and |x|>2. 
We then infer from the above bounds that, if CjW^p-'^^ > 1, 

1 + C7 iW'^f-^) ] _ 



Cs n(l,0) 

> / u{l, x) dx 

J{\x\>2} 

> \\u{l)\\i - / u{l,x) dx 

Jb(0.2) 



> M-[ [ \Vu{s,x)\'^ dxds -C8 2^ \\u{l)\\oo 

Jo Jrn 

/■I / \ N 

> M-Cs J [l + Cj M'^^P-^^ s") l|V'u(s)||«^ ds - Cg MP"^ 



> M-Cm M^^^f-^) M^"-?) s-w(A^+i) ds 







-10 
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the assumption q < p — {2N/(N + 1)) being used to obtain the last inequahty. Consequently, if 

C7M''(P-2) > 1, 

2^ Cg [C^ Af^^(P-2) - l] u(l, 0) > M - Cll (^M'liN(p-2)+2q) ^ ^pv^ _ 

Since prj G (0, 1) and r]{N{p — 2) + 2q) G (0, 1) by (p3|) we readily conclude that there are C12 and 
Ci3 such that ^(1,0) > Cu M'P'^ provided M > C13. 

We next consider an arbitrary initial condition uq fulfilling ([9]). Then there is C £ Cq°(]R^) such 
that < C < 1, supp C C B{0, 1), and 

C{x) uo{x) = ^ . 

5(0,1) ^ 

Denoting by u the solution to ([3|) with initial condition ■u(O) = ^ uq, it follows from the above 
analysis that u{l,0) > C12 (M/2)^^ if M > 2 C13. As the comparison principle warrants that 
^(1, 0) > 0), the expected result follows with C5 = 2 C13 and C4 = Cu 2"?"?. □ 



Proof of Theorem 14], For A > 0, t G [0, 00) and x G M , we define 

ux{t,x) := A^-"^ u {X^''~Pt,xo + X^^^+^x) , 

and observe that ux also solves ([3]) with initial condition ux{0). Furthermore, ux{0) fulfils Q and, 
if A > ^-VCp-i--?)^ we infer from ([22]) that 

[ ux{0,x)dx = XP-i+N{p-i~i) I uo{x)dx 

> A XP-'>+^^P-^-'>^ [ \x - xol'^P-'i^/^P-^-''^ dx 

> CuA. 

Therefore, if A > J-VCp-i-?) and A > Aq := C5/CU, we infer from Lemma [TOl that 

AP-« u{X^'^-P, xo) = ux{l, 0) > C4 ( [ ux{0, x) dx\ > . 

\JBio,i) J 

Consequently, u{t,xo) > and thus xq G for t G (o, . We finally use the time 

monotonicity of 'P{t) established in Proposition [5] to complete the proof of Theorem HI □ 

Combining Theorem|3]with a comparison argument allows us to extend Theorem|3]to q G [q2,p—l) 
under stronger assumptions on the initial data. 

Proposition 11 Consider xq G OVq and assume that there are ^ > 0, 5 > 0, and r G (1,92) such 
that 

uo{x) >A\x- xo|(P~^)/(P-^-^) for x G B{xo, 6) . (39) 
Then, if q £ [q2,P — have xq G V{t) for all t > 0. 



14 



Proof. We put Aq := l/||Vuo|loo and U{t,x) := Ag"'' u (^Aq^"^*, Ag~^+^a;j for (i,x) € [0, oo) x 

Then U solves ([3]) with initial condition Uq := C/(0) and ||VJ7o||oo = 1- In addition, Xq := Xq~^~'^xo 
clearly belongs to the boundary of the positivity set of Uq and it follows from (j39p that 

Uo{x) < A\x- Xo|(P-'')/(P-i-^') for x G S [xQ^dXl"^-") . (40) 

Introducing ri := (r + (72)/2 € (r, §2) we deduce from (|l0|) that, if (5i G ^0, ^Aq"^"*^^ and x G 
B{Xo,6i), 

Uo{x) > Xl^-^-"'^/^"-^-''^ A \x - Xo\^P-''^/^P-^-''^ |^_Xo|('^-'-i)/«P-^-'-)(P-i"'^i)) 
> ^(^-'•i)/((p-i-'')(p-i-''i)) ^ i^, _ Xol^P"''!)/^^"^"'"!) 

whence 

Uo{x) < ^1 |x - Xol^P-'-'^/^P-^-'-'^ for xeB (Xq, <5i) (41) 

with Ai := Aq"^ ^ ""^^^^ ^ ""^ (jj*" ^ ^ ''^^^ ^. Furthermore, we can choose 61 sufficiently 

small so that A > Aj.^, the constant being defined in Theorem [H 

Now, on the one hand, it follows from (j4ip and Theorem [5] that the solution u to 

dtu - ApU + \Vup = , {t,x) eQ , 
u{0) = Uo, X G , 

satisfies u{t,XQ) > for t > 0. On the other hand, since ||VC/o||cx) = 1) we infer from ([3]) and (fT2|l 
that 

dtU - ApU + iVC/j^'i > dtU - ApU + \\VU\\2,~'' iVf/l"? > dtU - ApU + IVJ/]" = 

in Q. The comparison principle [20l Theorem 8.2] then implies that U{t,x) > u{t,x) for (t, x) G 
[0,00) X . Consequently, 

u{t,xo) = Xl^ U [xl-'%Xo) > Xy u (xl-^\Xo) > 
for t > 0. □ 



6 Comments, extensions, and open problems 

Theorem [2] gives a precise description of the behaviour for large times of solutions to ([3]) for p > 2, 
g G (l,p — 1), and non-negative and compactly supported initial data. It then gives rise to several 
related questions, concerning the situation when either the parameters p and q or the initial data 
do not fulfil the above conditions. 

In particular, a natural question is whether a similar result is valid for non- negative initial data 
which decay rapidly at infinity but are not compactly supported. A first step in that direction 
would be to identify a class of initial data for which the solutions to ([3]), (jH) satisfy the temporal 
decay estimates (fT3]l . 
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Returning to compactly supported initial data, it is tempting to investigate what happens when q 
reaches the boundary of the range — 1) analyzed here. On the one hand, the critical exponent 
q = qi = p — 1 with p > 2 offers an interesting study of matched asymptotics that we plan to 
describe in a future publication. On the other hand, the critical case q = I seems to be quite open 
and interesting, even in the semilinear case [131 [14] . 

For larger values of q > p — 1 the situation is expected to be more classical: asymptotic simplifi- 
cation should take place for q > q^, = p — {N/{N + 1)) in the sense that the large time behaviour 
will be governed by the diffusive part of the equation. For intermediate values of g G {p—^, Q*), the 
existence of very singular solutions has been established in [33] and they are expected to describe 
the large time behaviour of compactly supported solutions. 

Concerning waiting times, an interesting technical open problem is to figure out whether Theo- 
rem H] is also true for g € [q2,Qi) {Q2 is defined in formula (j23p ). 

As a final comment, let us mention that the techniques of this paper could possibly be applied to 
similar equations involving diffusion and absorption like 

QiU = An™ - |Vnj« in Q, 

but we have noticed that there are a number of difficulties, maybe technical. 

Equations with variable coefficients are also worth considering, as well as more general equations 
of the form ([1]). In fact, the area of asymptotic behaviour of parabolic equations with variable coef- 
ficients is quite open even in the linear case. We refer to [24] for an interesting recent development 
for linear parabolic equations. 
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